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ABSTRACT 



^ ■ We consider a measure of the peculiar velocity field derived from the Mark III com- 



^ ■ pilation of 2900 spiral galaxies (Willick et al. 1996b), using an analysis method that is 
substantially free of bias (Nusser & Davis 1995). We expand the velocity field in a set of 
orthogonal, smooth modes, reducing the data to a set of 56 coefficients fitted to a maximum 
^ ^ redshift of 6000 kms~^, and maximum spherical harmonic of Z = 3. The radial resolution of 
^ ■ the modes degrades with redshift, from 800 kms~^ locally to 3000 kms~^ at 4000 kms~^ 
. redshift. Equivalent mode coefficients can be computed for the gravity field derived from 
any whole-sky redshift catalog of galaxies, such as the IRAS 1.2 Jy survey (Fisher et al. 
1995). Given the coefficients of the expansions, one can compare the velocity and gravity 
^ ' fields on a galaxy by galaxy basis, or on a mode by mode basis. Detailed comparison 
■ shows the two independent fields to be remarkably aligned in general. There are, however, 
systematic discrepancies in the fields that lead to considerable coherence in the residuals 
between them. These residuals take the form of a dipole field in the LG frame that grows 
with distance; it is not consistent with a bulk flow residual. 

We perform a likelihood analysis in the mode-mode comparison to determine which 
value of /9 = 0*^'^/6 for the IRAS gravity field is the best fit to the Mark III velocity field, 
considering the errors and covariance in both the velocity and gravity coefficients. We find 
that the most likely value lies in the range /3 = 0.4-0.6. However, in contrast with results 
we obtain using simulated galaxy catalogs, the per degree of freedom for the fit is well 
in excess of unity, primarily because of the coherent dipole residuals at cz ^ 3000 kms~^. 
Thus, despite the general alignment of the Mark III velocity and IRAS gravity fields, they 



* The manuscript with all its embedded figures (many in color) can be retrieved from 
ftp://magichean.herkeley.edu/puh/marc/m3itf /vitf.ps.g^ 
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do not agree in detail, precluding a firm determination of (3 from these data sets at present. 
The method is capable of measuring (3 to an accuracy of 10%, but without understanding 
these systematic discrepancies, we cannot infer a value of (3 from these data. 

Subject headings: cosmology — dark matter — galaxies: clustering — galaxies: formation 
— gravity — large-scale structure of the universe 
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1. INTRODUCTION 



The observed large scale flows of galaxies are a direct probe of the nature of the under- 
lying dark matter in the Universe. Under the assumption that the large scale structure has 
formed via gravitational amplification of small initial density fluctuations, the observed ve- 
locity field yields a prediction for the present matter-density fluctuations given an assumed 
value for O (Nusser et al. 1991). Furthermore, quasi-linear gravitational instability theory 
can be used to reconstruct the primordial density fiuctuations from the observed fiows. 
This enables us to address fundamental questions about the nature of the initial fluctua- 
tions and their relation to Q,. For example, Nusser & Dekel (1993) have demonstrated that 
unless Q. is around unity the initial density fluctuations were non-gaussian. 

The largest virialized systems in the Universe are the rare, rich clusters of galaxies, 
with typical sizes of less than 1 Mpc. The large scale flows, on the other hand, are 
motions with coherence length of order 10 - 40 Mpc, more than an order of magnitude 
larger than the rich cluster sizes. Measurements of mass by means of the alignment which 
is expected between velocity and gravity flelds on these large scales is thus quite distinct 
from the virial, or hydrostatic, or weak lensing analyses that can be applied to the clusters. 
A trend of increasing mass to light ratio is well documented to exist with increasing scale 
size of the measurement (Ostriker et al. 1974, Davis et al. 1980). Rich clusters fall on 
this trend, with an inferred density parameter of order ?a 0.1 — 0.2, consistent with the 
baryon fraction in the clusters (e.g. White et al. 1993). A fundamental question is whether 
the trend of versus scale asymptotes at this value, or continues to climb upward on 
larger scales. If galaxy formation is more efficient in denser regions or if galaxies have 
some degree of velocity bias, or if the dark matter is too hot to fully fall into the cluster 
potential wells, then we might expect to measure a larger value of VLq on the scale of the 
large-scale flows. Considerations of an inflationary early Universe have traditionally led 
to the expectation of ilo = 1 (although see e.g. Liddle and Lyth 1993, Bucher and Turok 
1995); our confidence in such speculations would be considerably enhanced if someone 
could actually produce compelling evidence, on any scale, to support such a high density. 

On large scales, virial equilibrium does not apply, but one can still derive estimates of 
the mass density, now based on fiow models, either linear or nonlinear. The earliest detailed 
fiow models were based on Virgocentric infall (e.g. Aaronson et al. 1982, Tonry and Davis 
1981). These models assumed a smooth background and spherical symmetry for the local 
super cluster, and this restrictive symmetry allowed the incorporation of nonlinearity in 
the relationship between velocity and density. We now know that spherical symmetry 
is a rather poor approximation for the infall to any cluster, but linear theory analysis, 
which assumes no special symmetry, remains a valid tool. The large scale fiows are almost 
certainly the result of the process of gravitational instability (e.g. Dekel et al. 1993, but see 
Babul et al. 1994), with the overdense regions attracting material, and underdense regions 
repelling material. Initial conditions in the early universe might have been somewhat 
chaotic, so that the original peculiar velocity field (i.e. deviations from Hubble fiow) was 
uncorrelated with the mass distribution, or even contained vorticity. But those components 
of the velocity field which are not coherent with the density fluctuations will adiabatically 
decay as the Universe expands, and so at late times one expects the velocity fleld to be 
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aligned with the gravity field, at least in the limit of small amplitude fluctuations. In this 
paper we shall attempt to make this comparison directly, using the best presently available 
data for both the velocity and gravity fields. 

The subject of large scale flows has been extremely active of late, and two excellent 
reviews are now available (Dekel 1994; Strauss and Willick 1995). The observational 
situation is rapidly improving; in this paper we shall consider the velocity held derived 
from a sample of 2900 spiral galaxies with TuUy- Fisher measurements (Willick et al. 1995; 
Willick et al. 1996a,b). The basic principle is to use late-time linear theory to relate 
the measured velocity fleld to the gravity fleld as predicted from the observed galaxy 
distribution. Given complete knowledge of the mass fluctuation fleld 6p{r) over all space, 
one could trivially write the fluctuations in the gravity fleld g(r) as 



where p is the mean mass density of the Universe. The great simpliflcation of late time 
linear theory is the intuitive relation between the gravity fleld g and the velocity fleld vl, 
namely vl = g t where the only possible time t is the Hubble time. The correct expression 
(Peebles 1980) is 



If the galaxy distribution at least approximately traces the mass on large scale, with linear 
bias b between the galaxy fluctuations Sq and the mass fluctuations (i.e. Sq — bSp), then 
from (1) and (2) we have 



where n is the true mean galaxy density in the sample, /? = and where we have 

replaced the integral over space with a sum over the galaxies in a catalog, with radial 
selection function (p{r) . {(p{r) is deflned as the fraction of the luminosity distribution 
function observable at distance r for a given flux limit; see e.g. (Yahil et al. 1991).) The 
second term is the correction for the uniform component of the galaxy distribution which 
would exactly cancel the flrst term in the absence of clustering within the survey volume. 
One problem with Equation (3) is that the sum is to be computed in real space, whereas 
the galaxy catalog exists in redshift space. Furthermore, the integral is to be computed 
over all space, but the sum over the catalog must be truncated at some flnite radius. Note 
that the result is insensitive to the value of Hq, as the right hand side has units of velocity. 
We shall henceforth quote all distances in units of kms~^. The density n and selection 
function 0(r) can be estimated from the redshift catalog used in the analysis, so that 
comparison of the measured velocity fleld to the predicted velocity fleld VL(r) gives us a 
measure of /3. The degeneracy between O and b can be removed in principle by nonlinear 
effects (Nusser et al. 1991, Dekel et al. 1993), but in practice one should not be optimistic 
because, among other problems, nonlinear bias will enter at the same order as nonlinear 
velocity effects. 




(1) 




(3) 
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The direct comparison of the pecuhar velocity of selected galaxies versus the gravity 
predicted from Equation (3) is fraught with difficulty. Distances to individual galaxies 
arc typically uncertain at the level of 15-20% and are furthermore subject to considerable 
Malmquist bias. To compute the gravity field, one needs a nearly whole sky redshift catalog 
of galaxies. To date the best available catalog is the 1.2 Jy flux limited IRAS catalog of 
5300 galaxies (Fisher et al. 1995), which we shall use in this analysis. The IRAS predicted 
gravity field has substantial uncertainty because of the dilute sampling of the density 
field. Furthermore, the predicted velocities are subject to large covariance, because the 
IRAS predicted field is quite non-local. Flows within the turnaround radius of large clusters 
are multivalued, so that in projection a given redshift can correspond to three separate 
distances. The examination of numerical simulations demonstrates that coherent shell 
crossing is a common occurrence near dense centers: the redshifts of a set of points can 
be larger (smaller) than the redshifts of nearby points with smaller (larger) distances. In 
this case, essential information has been lost and recovery of the correct distances from 
redshift information alone is impossible. 

Several comparisons of velocity and gravity fields have been made with older datasets 
(e.g. Hudson 1994; Kaiser et al. 1991; Yahil 1988; Strauss and Davis 1988) but these 
analyses were all meant to be preliminary and none of their conclusions were compelling. 
To date nobody has included a proper treatment of the correlated noise in the analysis. 
Shaya et al. (1995) have compared the peculiar velocities of a set of 298 spiral galaxies 
within a redshift of 3000 kms~^ to the gravity field derived by least action analysis of 
the 1138 mass tracers dervived from the Nearby Galaxies Catalog (TuUy 1988). Their 
preliminary analysis gives a value fio = 0.2 ± 0.2, but they have not included covariance 
in the errors of the predicted mass model, nor have they realistically included possible 
infiuence from the mass distribution at redshifts z > 3000 kms~^, nor is it clear how well 
their analysis compares with the linear theory techniques used by others. 

Because of the covariance of the IRAS gravity field a comparison of the mass fluc- 
tuations as inferred from the flows versus the IRAS density fleld has considerable merits. 
Assuming the flow is irrotational, the POTENT algorithm (Dekel et al. 1990) reconstructs 
a three dimensional flow from the observed radial flow field, and then computes derivatives 
of this fiow to infer the underlying mass field (Nusser et al. 1991). This then can be com- 
pared directly to the IRAS density field (or to any other survey). Such a comparison was 
done by Dekel et al. (1993), using an earlier version of the IRAS catalog, and the Mark II 
catalog of peculiar velocities for 493 objects (Faber and Burstein 1988). On the basis of 
this density- density comparison, Dekel et al. (1990) derived (3 = 1.3 ± 0.3. More recently, 
Hudson et al. (1995) have compared POTENT mass densities reconstructed from the 
Mark III data (Willick et al. 1995, 1996a, 1996b) with the optical density field of Hudson 
(1994), finding /3 = 0.74 ±0.13. A new analysis using the 1.2 Jy IRAS survey and the 
Mark III peculiar velocity catalog is in preparation (Dekel et al. 1996). 

Given the significance of the questions raised, all viable alternative algorithms should 
be pursued to extract the implications of the large scale flows. In this paper we shall 
compare the observed radial velocity field with the radial gravity field inferred from the 
1.2 Jy IRAS redshift survey (Fisher et al. 1995). We shall use the method of orthogonal 
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mode expansion as described by Nusser and Davis (1995), which we shall henceforth label as 
the Inverse Tully-Fisher (ITF) method. A disadvantage of the density-density comparison 
is that the mass-fluctuations which are derived by POTENT increase with decreasing p. 
According to the quasi-linear approximation (Nusser et al. 1991) employed by POTENT, 
the mass-fluctuation fleld inferred from the velocities diverges for low thus invalidating 
the density-density comparison. Therefore the density-density comparison can not formally 
rule out very low values of 13. However, the fluctuations in the velocity fleld predicted from 
the observed distribution of galaxies are smaller for lower /3 and the applicability of the 
linear approximation is guaranteed. Therefore, the velocity-velocity comparison provides 
a better tool for assessing the agreement between the observations for low (3. 

This paper closely follows the methods described in Nusser and Davis (1994, hereafter 
ND94) and in Nusser and Davis (1995) (hereafter ND95), with some additional refinements 
developed since then (Davis and Nusser 1995). In Section 2, we review the essentials of 
these methods and describe our choice of basis functions. We also discuss the procedure 
for fitting the modes to both the velocity and gravity fields. As has been the case in the 
past, mock catalogs constructed from A^-body simulations are essential for debugging and 
calibrating the methods. This is especially so for our application, since the entire analysis 
is performed in essentially pure redshift space. In Section 3, we show a few new details 
on the tests applied to the mock catalogs. Section 4 gives details of the comparison of the 
Mark III velocity field with the IRAS inferred gravity field, on a point by point basis, as 
well as on a mode by mode basis. A detailed likelihood analysis is presented in Section 5, 
where the needed covariance matrices are derived as well. Section 6 presents a discussion 
and final conclusions are given in Section 7. 

2. Reconstruction of Peculiar Velocities 

In this section we outline the methods of ND94 and ND95 for deriving the smooth 
peculiar velocities of galaxies from an observed distribution of galaxies in redshift space 
and, independently, from a sample of spiral galaxies with measured circular velocities 
Tj = logAt) — 2.5 and apparent magnitudes rn. Any of the available methods (e.g. Yahil et 
al. 1991, Fisher et al. 1995b) for generating IRAS predicted peculiar velocities can be used. 
However, the method of ND94 is particularly convenient, as it is easy to implement, fast, 
and requires no iterations. Most importantly, this redshift space analysis closely parallels 
the ITF estimate described below. Since we shall use this analysis in what follows, we next 
present a very brief summary of the method. 

2. 1 Peculiar Velocities from the Distribution of Galaxies in Redshift Space 

Let v{s) be the comoving peculiar velocity relative to the motion of the Local Group 
(LG) and s the comoving redshift space coordinate. (We use s = cz/Hq, noting that higher 
order terms in z have a negligible influence on the derived velocity fleld). To first order, 
the peculiar velocity is irrotational in redshift space and it can be derived from a potential 
function, $(s), defined by v{s) — — V$(s). Let 6o{s) be the fractional fiuctuation of the 
density field traced by the distribution of galaxies in redshift space. If we expand the 
angular dependence of $ and Sq in terms of spherical harmonics (Fisher, Scharf & Lahav 
1994, Fisher et al. 1995b) in the form, $(s) = X)^o T,L=-l ^im{s)Yim{d, (p) and similarly 
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for Sq, then, to first order, and Sqi^ satisfy. 



s'^ ds \ ds 



m 



l + f3 



1 Ki + i)^i 



m 




s dins ds 



1 dln0 d$/. 



m 



(4) 



where (/> is the selection function of the sample. As emphasized by ND94, the solutions 
to (4) for the monopole (/ = 0) and the dipolc (/ = 1) components of the radial peculiar 
velocity in the LG frame are uniquely determined by specifying vanishing velocity at the 
origin. That is, the radial velocity field at redshift s, when expanded to harmonic Z < 1, is 
not influenced by material at redshifts greater than s. 

As an estimate of the field 5q{s) from the discrete distribution of galaxies we consider. 



We emphasize that the coordinates s are in observed redshift space, expanded in a galactic 
reference frame. The only corrections from pure redshift space coordinates is the collapse 
of the fingers of god of the known rich clusters prior to the redshift space smoothing 
(Yahil et al. 1991). The density field is here smoothed by a gaussian window with a 
redshift dependent width, a"(s), which is proportional to the mean particle separation, i.e, 
a{s) oc (n(/)(s))~^/"^, where n is the mean density of the IRAS sample. The purpose of the 
smoothing is to limit the discreteness noise of the gravity field constructed from the dilute 
redshift sample. This particular smoothing ensures that the error in the field (5) due to 
Poisson noise is constant in space for the case of a uniform distribution of galaxies. A more 
elaborate realization of the density field based on Wiener filtering is described by Fisher 
et al. (1995b), with results very close to those obtained here. Weighting the galaxies in 
(5) by the selection function evaluated at their redshifts rather than the distances which 
are unknown yields a biased estimate for the density field and gives rise to Kaiser's rocket 
effect (Kaiser 1987). However, this bias should be compensated by the term involving the 
selection function in (4). 



Given a sample of galaxies with measured circular velocity parameters, iji, apparent 
magnitudes, rrii, and redshifts, Zi, the goal is to derive an estimate for the smooth un- 
derlying peculiar velocity field. We assume that the circular velocity parameter, 77, of a 
galaxy is, up to a random scatter, related to its absolute magnitude, M, by means of a 
linear inverse TuUy-Fisher (ITF) relation, i.e.. 



The advantages of inverse TF methods were first described in detail by Schechter 
(1980) and Aaronson et al. (1982). Most critically, samples selected by magnitude, as most 





2.2 Peculiar Velocities from the Inverse Tully-Fisher relation: 
Expansion of the Fields in Orthogonal Functions 



T] = 7M + r]o. 



(6) 



7 



are, will not be plagued by Malmquist bias effects when analyzed in the inverse direction. 
Methods based on an inverse relation generally require a velocity model (Schechter 1980, 
Aaronson et al. 1982, Lynden-Bell 1991, Lahav 1991) which can be parameterized in terms 
of observable redshifts. ND95 consider a very general model and write the absolute magni- 
tude of a galaxy, Mj = Mq^ -t-Pj, where Mqi = mi + b\og{zi) — 15 and Pi = 5log{l — Ui / Zi) , 
where rrii is the apparent magnitude of the galaxy, Zi is its redshift in units of kms~^, 
and Ui its radial peculiar velocity in the LG frame. In general, one can write the function 
Pi in terms of an expansion over orthogonal functions. 



J max 



J=0 



with orthonormality conditions. 



EW = ^^': (8) 

i=l 

with the zeroth mode defined by F° = l/\/iVg, where Ng is the number of galaxies in the 
sample. The zeroth mode describes a Hubble-like flow in the space of the data set which is 
clearly degenerate with the zero point of the ITF relation. Here we arbitrarily set = 0. 
The best fit parameters, , the slope, 7, and the zero point 770, are found by minimizing 
the statistic 



X 



,2 ^ ^ hMpi + 'yPi + r]o-r]i) 



where o"^ is the rms scatter in ry about the ITF relation. With our choice for the zeroth 
mode, the solution for is. 



= --J2{lMoi-Vi)F^, 

^ i 



Em' ^OiMo,. (1 + Ng E,- FiFl) - Ng M, 



As described in ND95 the choice of radial basis functions for the expansion of the 
modes can be made with considerable latitude. The functions should obviously be lin- 
early independent, and close to orthogonal when integrated over volume. They should be 
smooth and close to a complete set of functions up to a given resolution limit. ND95 chose 
spherical harmonics Yj^ for the angular wavefunctions and the derivatives of spherical 
Bessel functions for the radial basis functions, motivated by the desire to use functions 
which automatically satisfy potential theory boundary conditions at the origin and the 
outer boundary. That is, they chose 



m=l 



P{Z, ^» = E E E ^ (^'l ^^nZ) - Cn) 0) . (11) 



2; 

n=0 /=0 m=—l 
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The constant cn is non-zero for the dipole term only and ensures that P = at the origin, 
and is non-zero at the outer boundary. Details of how the orthogonalized functions F- are 
derived from this expansion are given in ND95. 

To keep the expansion of the basis functions compact, it is sensible to match the 
resolution of the wavefunctions to the spatial distribution of the TuUy-Fisher data. Since 
the precision of the inferred peculiar velocities is a constant fraction of the distance to any 
galaxy, spherical Bessel functions, which oscillate uniformly with their argument, are not 
the optimal choice of wavefunction. Furthermore, the mean intcrparticle spacing of galaxies 
in flux limited redshift catalogs, from which we infer the density and gravity fields, can 
increase even more rapidly than linearly with distance. For example, for z ^ 6000 kms~^, 
the mean intcrparticle spacing within the 1.2 Jy IRAS survey increases as z^'^. Therefore, 
to make optimal use of the available information from the ITF field, we expand the modes 
of Pi in terms of a transformed radial coordinate that decreases the oscillation periods of 
the Bessel functions roughly according to the distribution of data as a function of redshift. 
A convenient, but by no means unique choice is 

y={ln[l + {z/z.r]f' . (12) 

For an increment Ay, this function has resolution Az = {dy/dz)~^Ay that is constant 
on small scale and Az oc zln{z) on large scale. We set z^ = 1000 kms~^. The radial 
wavenumbers for each / value have been chosen for the lowest values of kn such that 

jl-liknVmax) = . (13) 

This choice does not restrict the value of the expanded velocity field near the boundary 
Umax- As we shall see, transforming to the variable y results in a significant reduction in 
the value of the of the fit relative to the expansion in terms z. Given these choices 
of radial variable and values of k^, the resulting wavefunctions for / = 0,4 are plotted in 
Figure 1. Note how the wavelength of the modes increases with distance. All of the modes 
are forced to be orthogonal to the P — constant mode when summed over the set of points 
in the Mark III catalog, and this is accomplished by small additive constants. However, 
the n = 1, I = mode is positive definite over the volume, and to make it orthogonal 
to the P = constant mode would require a substantial additive offset, which would then 
violate the proper boundary condition at the origin. We therefore delete this mode from 
inclusion in our expansion set. 

As a demonstration of the use of the variable y, we tested the Mark III catalog which 
lists 2237 galaxies within s — 6000 kms~^, details of which will be described below. We 
took the linewidth scatter to be cr^ = 0.05 (Willick et al. 1995). For no velocity model, 
(i.e. setting all coefficients = 0), the scatter diagram of the TuUy-Fishcr regression 
of observed versus predicted linewidths r] yields = 3516. Fitting a flow model of 69 
degrees of freedom {Imax = 3, nmax = 4) in terms of radial wavefunctions i[{knz)/z leads 
to a much reduced — 2667, while doing the same fit of 69 modes in terms of ji{kny)/z 
gives = 2644, a modest but significant improvement. Further tests on simulations with 
smaller cr^ demonstrate that the effect of the coordinate transformation is considerably 
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more pronounced. There is clearly room for refinement here, but this choice of radial 
variable is sufficient for our purposes, especially since we shall use identical basis functions 
for the expansion of the velocity and gravity fields. 

2.2.1 Choosing the resolution limit 

In the comparison of the Mark III velocity and IRAS gravity fields discussed below, we 
shall use the Mark III sample only within a redshift limit of 6000 kms~^. This corresponds 
to setting the outer boundary condition at Umax = 1-9. If we consider that the number of 
peaks in the nth radial wavefunction is n, then the effective resolution R at any redshift 

is R{z) ^ (^g) For n^ax = 4, Zmax = 6000 kms-^ and z^ = 1000 kms-^ the 

resulting limiting radial resolution is 680 kms~^ at z = 1000 kms~^, and 3400 kms~^ at 
z = 4000 kms-^ 

The angular resolution for a spherical harmonic function Y/^ is 180// degrees. For 
redshift z, this corresponds to a maximum spatial resolution of i? ~ Trz/lmax- If we desire 
to have the angular and physical resolutions approximately equal at 3000 kms~^, where 
the radial resolution is approximately 2400 kms~^, then we should select Imax = 3.9 ~ 4. 
We shall in fact limit Imax to 3, in order to reduce the number of degrees of freedom. The 
result is that our radial resolution is slightly higher than our transverse resolution, but by 
our choice of radial variable y, the ratio of transverse to radial resolution is approximately 
constant with distance, which is a desirable feature. 

The choice of the variable y is governed by the desire to have a mode resolution that 
is a fixed fraction of the resolution of the TF data, but as long as the resolution is below 
that of the intrinsic data, we can be assured that both the smoothed velocity and gravity 
fields have identical resolution. We have chosen a mode resolution corresponding to the 
spatial frequency of the fourth dipole mode: k = That is, we include the first 4 radial 
modes that satisfy the outer boundary condition, equation (13), for the harmonic 1 = 1. 
We use this limiting cut in wavenumber as a limit for the other spherical harmonics. Thus 
to the same limiting value of /c, we should include 4 modes for / = 0, the first of which is 
discarded, 3 modes for / = 2, 3 modes for / = 3, 2 modes for / = 4, 2 modes for / = 5, and 1 
mode for / = 6. However, tests based on the efficiency of the decrement in the value of 
the fit demonstrated that a cutoff in I is desirable here. Thus on top of the sharp cutoff in 
/c-space we eliminate modes with harmonic index exceeding Imax = 3, as mentioned above. 
This yields 3x1 + 4x3 + 3x5 + 3x7 + 5 = 56 modes in total. The Mark III velocities 
within a limiting redshift of 6000 kms~^ were then derived from the ITF method. The 
peculiar velocities as predicted from the IRAS real or mock catalogs were expanded in the 
same ITF basis functions, as described next. 

2.2.2 Expanding the IRAS velocity field 

Using the machinery for computing a gravity field described in §2.1, one can generate 
a linear theory predicted peculiar velocity vl for any point in space as a function of 
its redshift for any value of We must ensure that the smoothing scales of the ITF 
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and IRAS predicted peculiar velocities are matched to the same resolution. Therefore 
we expand in terms of the modes used in the velocity model (12). Because of the 
orthonormality, we can write the mode coefficients as 



where the H' term is a correction for Hubble flow and the summation index i is restricted 
to be over the positions of the same galaxies in the ITF expansion. 

This procedure will fllter out higher frequency modes that are not described by the 
resolution of our basis functions. However, it is important to recall that the modes are 
chosen to be orthogonal to the mode Pi = constant, which would describe a smooth Hubble 
flow. In the fltting for the ITF modes, pure Hubble flow is absorbed into a shift of the 
zero point ?7o ^-nd the orthogonality is ensured. Within a given set of test points occupying 
a volume smaller than that used to deflne the gravity fleld, it is possible for to have a 
non-zero value of Hubble flow H' , which must be removed from vl before we tabulate the 
mode coefficients. That is, we tabulate the mean Hubble ratio 



and subtract it from the predicted field vl- This "breathing mode" which mimics a Hubble 
fiow is not trivial in amplitude, and can be a 10% correction on the effective Hubble 
constant within simulated catalogs. 



ND94 and ND95 tested the velocity reconstruction methods and their biases using 
mock catalogs drawn from an A'"-body simulation of a CDM, {h = 0.5, flo = 1, and erg = 
0.7) universe. The simulation had the disadvantage of having low amplitude large scale 
flows and therefore was not fully suitable for our purposes. Furthermore, the algorithm 
for attenuating the rather large small-scale incoherent velocities seen in the simulation 
resulted in too much reduction of their amplitude. Here we provide more powerful tests 
of the methods using mock Mark III and 1.2Jy IRAS catalogs which are extracted from 
a pure n = —1, Oq = 1 simulation of 10^ points kindly supplied by S. White. We use 
a normalization as = 0.7, which then determines that the size of the box is L = 235h~^ 
Mpc. Such a simulation has the desirable feature of having very high-amplitude large scale 
flows. In this simulation, as in the CDM simulation used in ND, the one- dimensional small 
scale (few Mpc scale) velocities are far in excess of that observed for the real universe 
(fa 140 — 200 kms~^; see Miller et al. 1996) and have to be smoothed somehow. Instead 
of replacing the A^-body peculiar velocities with Gaussian smoothed velocities as we have 
done in previous work (Fisher et al. 1993), we replace the original A^-body velocities with 
the average of a Gaussian smoothed fleld (with smoothing length a = 100 kms~^) and the 
original field. Thus we reduce the one dimensionsal RMS amplitude of the very small scale 
velocity field cr^ from 400 to 200 kms~^ while keeping the large scale velocities unchanged. 




(14) 




(15) 



3. Tests with Mock Catalogs 
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This is of course not a self-consistent procedure, but better approximates the observations 
and the expected effects of possible velocity bias, or of baryon dissipation within extended 
dark matter halos. In dense regions, with 6 > 600, we collapse the redshift space fingers 
of god, as is done in the centers of real clusters for both the Mark III and IRAS catalogs. 

In order to generate realistic mock catalogs from the simulation, a particle moving at 
600 kms~^, with small local shear, and which is not in a dense region is chosen as the 
central observer. A mock full sky 1.2Jy IRAS catalog is generated from the simulation 
by assigning luminosities to particles according to the luminosity function of the 1.2Jy 
IRAS survey. This procedure results in 5021 "galaxies" out to a redshift of 180h~-'^Mpc 
from the "observer" . 

To select mock Mark III catalogs, we follow the procedure described in ND95, with 
the following modifications: we extract a nonuniform subset of points from the mock 
IRAS catalogs, with galaxies at |6| < 18° heavily selected against, and with objects in 
the hemisphere of the microwave dipole anisotropy vector having twice the probability of 
selection as objects in the anticenter direction. Each point within the mock catalog has 
an absolute magnitude drawn from the luminosity function of the 1.2 Jy IRAS catalog. 
The selected points in the mock Mark III sample are given a linewidth i] drawn from a 
linear TuUy-Fisher law with ryo = —2 and 7 = —0.1, and with gaussian scatter — 0.05, 
a value we find to closely fit the observed scatter in the real Mark III catalog. The mock 
Mark III data thus consist of a position on the sky, a redshift, apparent magnitude, and 
linewidth for each galaxy; the true peculiar velocity of each "galaxy" is used only to test 
the reliability of the reconstruction. A much more elaborate construction of mock Mark 
III catalogs is described by Kolatt et al. (1996). We have tested our procedures on this 
catalog for which we have no more difficulty recovering the true underlying field than with 
the mock catalogs described here. 

With these two catalogs, we can generate mock ITF velocity fields as described above, 
and mock IRAS velocity fields using the technique of ND94 described above. By simply 
changing the random number seed, we have generated different realizations of the mock 
Mark III sample for the same points, and different mock IRAS catalogs from the same 
central observer within the A'^-body simulation. 

To generate the mock IRAS velocity field by means of the method described in Section 
2.1, we first use the unbiased "galaxy" distribution out to 18000 kms~^ to generate a 
smooth density field. We fix the Gaussian smoothing length a{s) to be the mean particle 
separation at the redshift s. For the 1.2 Jy IRAS catalog, this smoothing is approximately 
cr(s) = 345 kms~^ at ^ = 1000 kms~^, a{s) = 570 kms~^ at 2 = 3000 kms~^, and 
(t{s) = 980 kms~^ aX z — 6000 kms~^, considerably smaller than the mode filtering 
we shall apply to the ITF velocity and IRAS gravity fields, but large enough to filter 
out many of the nonlinear structures in the mock IRAS catalog. The density field and 
resulting velocity field (equations 4 and 5) are then expanded in spherical harmonics up 
to Imax = 12, although in fact a lower cutoff at / = 4 is sufficient. Using this relatively 
high resolution velocity field based on the mock IRAS catalog, we generate a linear theory 
prediction for each point in the mock Mark III catalog. 
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We use the 56 mode expansion to minimize the scatter of the TuUy-Fisher relation 

in the mock Mark III catalog; for the 1733 points in the mock catalog within s = 6000 
km/s, the of the fit is reduced from 2264 with no fiow, to 1900 with a 56 mode fiow 
model. Note that the reduced is somewhat greater than unity in spite of our precise 
knowledge of the scatter in 77. The slightly high is simply due to components of the 
velocity field with higher spatial frequency than the resolution of our mode expansion; 
however, increasing the number of modes in the expansion reduces x^ relatively slowly. If 
the r] scatter were smaller, more degrees of freedom would be worthwhile. We recover best 
values of r]o = —1.94 and 7 = —0.0965, only slightly off their input values of —2.0, —0.1. 

We find it most useful to present the results of the fiow model pictorially. In Figure 2 
we present three redshift slices of the ITF velocity field unf derived for the mock catalog 
using the 56 mode expansion. The plots are in galactic coordinates where stars (circles) 
represent mock galaxies outfiowing (inflowing) in the LG frame, with sizes proportional 
to the fiow amplitude. In Figure 3 we show uuf projected onto the Supergalactic Plane. 
Note that although the data are irregularly distributed, the velocity field is very smooth, 
by construction. These plots, like all the plots in this paper, are in the Local Group frame. 
Note the strong reflex dipole signature which dominates the large scale flow. 

The real advantage of working with mock catalogs derived from self-consistent A^-body 
simulations is that for each point we know the true peculiar velocity Utmei including all the 
nonlinear, small scale effects. Figure 4a shows the true, unfiltered velocity field Wtrue of the 
mock Mark III points, projected to the supergalactic plane, while Figure 4b shows the true 
field after filtering by the ITF mode functions Wtrwe' Equation (14). Note how the region 
of substantial shell crossing shown near the clusters of the upper left section of Figure 4a 
is largely eliminated in the filtered version of the field, Figure 4b. Figure 4c shows the 
IRAS derived, unfiltered velocities Uiras-, while Figure 4d shows the IRAS derived velocities 
after filtering, u{^^g. The Utme field has a dipole signature on large scale, and considerable 
complex structure around the dense knots; the Uiras field is considerably smoother and 
cleanly separates the infalling and outflowing zones, particularly in the upper left portion 
of Figure 4c, as expected from a Zeldovich approximation methodology. Note that Figures 
3, 4b, and 4d, the smoothed versions of the ITF, true, and IRAS fields, are all very similar, 
indicating that the ITF process does recover the true peculiar velocity field, as does the 
IRAS gravity field. The similarity of Figure 4c to 4d demonstrates that the 56 mode 
filtering recovers most of the information contained within the IRAS predicted flow field. 
The higher resolution velocity information of the unfiltered IRAS flow fleld has relatively 
low amplitude in most regions. 

To be slightly more quantitative, we present in Figure 5a the residual unf — wfrwe ^ 
a function of z, with points within 60 degrees of the apex angle toward (/ = 270, b = 30) 
plotted on the right, and those within 60 degrees of the antiapex plotted on the left. Note 
that there is a slight monopole residual, with all distant points having slightly negative 
residuals. This has resulted from bias in the determination of the TuUy-Fisher slope 7, 
which has been fit by the x^ minimization procedure but has led to a biased estimate, 
7 = —0.0965, instead of the correct value 7 = —0.1. The bias has been induced by the 
thermally "warm" velocity field of the mock Mark III catalog, as expected (see ND95 for 
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further details). If instead of allowing 7 to be determined by the minimization, we fix 
it at its correct value, the resulting increases from 1900 to 1905, but the resulting field 
has no residual monopole, as shown in Figure 5b. Figures 2-4, as shown, were corrected 
for this bias. We can expect a similar bias for the real data, since the scatter in Ou has 
approximately the same amplitude as the mock catalog, and we shall therefore correct the 
measured TF slope by the same factor of 3.5%. This bias is so minor that none of the 
results to follow are affected if we ignore it altogether. 

Having corrected for this bias, we next examine maps of the residual fields. Figure 
5c shows the residuals unf — plotted in the same manner as Figure 5ab. Note that 

the residuals are now somewhat larger and increase with distance, but show no systematic 
offsets. Figures 6a and 6b show the residuals u^f —ui^^^ and u^f — ui^g^g plotted in slices of 
the full sky. Note that the amplitude of Figure 6b is larger than in Figure 6a, showing that 
IRAS induced errors are a major contributor to the errors in the procedure. Note again that 
all the errors are coherent, as expected since the field is limited to / = 3; neither the velocity 
nor gravity field contain any information at higher resolution. Comparison of different 
mock IRAS and mock Mark III catalogs shows similar coherence in the residuals, but 
largely differing sky patterns than shown in Figure 6b. An important point to emphasize 
for what follows is that the dipolc component of the error in the mock catalogs is small; 
the errors appear to be largely dominated by higher multipoles. The coherence of the 
residuals makes it difficult to judge the adequacy of the fit, but we shall see below when 
we consider the statistics of the mode amplitudes that these fits are consistent with the 
expected noise. 

4. Comparison of Mark III vs IRAS 

Having set the stage by means of the mock catalogs, we now consider the actual 
comparison of the ITF velocity field from Mark III and the gravity field derived from the 
IRAS 1.2 Jy catalog. Although the agreement between the mock Mark III and the mock 
IRAS catalogs was outstanding, we shall find that the comparison of the true velocity and 
gravity fields is considerably more problematic. 

4.1 Preparation of the Mark III Database 

Before constructing the ITF velocity field, it is necessary to bring the Mark III data 
to a uniform linewidth and magnitude system. In its raw form, the catalog consists of 
six distinct pieces, each with its own TF relation (Willick et al. 1996a), but by means of 
direct comparison of overlapping galaxies, the catalog can be standardized. We brieffy 
summarize the procedure here; complete details are given by Willick et al. (1996b). Note 
that we use only spiral galaxies in the present analysis. 

We assume that all six Mark III samples can be placed on a "common system" by 
linearly transforming the apparent magnitudes and and velocity widths, or TF observables. 
For the reasons discussed by Willick et al. (1995), the Han-Mould (HM) cluster sample 
was used as the template data set, i.e., the HM data are unchanged in our analysis. The 
remaining samples are then transformed via a hierarchical overlap-comparison procedure 
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analogous to that described by Willick et al. (1996a), except that now the TF observables 
rather than the TF distances are required to agree in the mean with those of the same 
objects appearing in aheady-transformed samples. Transformations of the form 

^/common = ^0,5 + Cli^sVS + 0'2,S'n% (16) 

and 

^common = "^5 + 60, S + ^l,s(w5 - 5 log C^) (17) 

are fitted to the overlap data (with one exception, discussed in the next paragraph) , where 
the subscript "S"' denotes a given sample, and "common" refers to the values of the 
observables on the "common system" ultimately adopted for our analysis. The quadratic 
coefficient 02,5 was non-zero only in the case of the Courteau-Faber (CF) sample, while 
the coefficient 61,5, which represents a luminosity dependence of galaxy color, is non-zero 
only in the case of the if-band Aar onson et al. (A82) sample. Values of all coefficients 
that appear in equations (15) and (16) are tabulated by Willick et al. (1996b). 

For the Mathewson et al. (1992; MAT) sample, we adopted a modification to Equa- 
tion (15). When the MAT ry's were initially transformed to the common system, the fit 
residuals correlated weakly with redshift, in the sense that higher redshift MAT objects 
appeared to have smaller widths than those of the same objects in other samples. A possi- 
ble interpretation is that the MAT widths, which were measured chiefiy in the 21 cm line 
at lower redshifts, but primarily optically in the Ha line at higher redshifts, are biased low 
with increasing redshift. To account for this possibility, we tried a transformation of the 
form 

?7common = ao,MAT + ai,MAT'75 + OS.MAT log(c2;). (18) 

This yielded a coefficient 03, mat = 0.065 ± 0.025, which is small but statistically signif- 
icant, and a small but significant decrease in the scatter of the overlap comparison. We 
subsequently found that the ITF velocity fit using the MAT widths transformed according 
to (17) produced a significant reduction in the ITF statistic. We have therefore adopted 
the redshift-dependent transformation (17) for the MAT sample throughout the analysis 
presented in this paper; we discuss possible implications of this further in § 6. 

A final step in the preparation of our modified Mark III sample involves the rejection 
of outliers. As has been noted elsewhere (e.g., Willick et al. 1996a), the TF residuals 
of typical complete samples are not strictly Gaussian. Rather, they tend to consist of a 
gaussian core plus a small (1-2%) admixture of 4-5(t outliers. Such objects are not helpful 
to velocity analyses, as they vitiate the assumption of Gaussianity that underlies almost 
all statistical methodologies. In our study we purged such points in the following way. 
First, we fit a bulk flow model to the entire sample. We computed the scatter cr^ from 
this fit, and then eliminated the 29 objects which deviated by more than 3.3a"^. These 
rejected galaxies were distributed throughout the volume of the MAT sample. When we 
refitted the bulk fiow models, the scatter was reduced by 10% as a result of this small 
pruning of the sample, while the direction and amplitude of the bulk fiow were essentially 
unchanged. Moreover, there were no remaining points which deviated by > S.Scr^j relative 
to the reduced scatter. Thus, our exclusion procedure resulted in a statistically improved 
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(i.e., more gaussian) sample, without significantly altering the global properties of the 
velocity field the sample produces. 

4-2 Maps and Residuals of the Real Sky 

The first stage is to fit the TuUy- Fisher data of Mark III to the 56 orthogonal modes. 
For 2237 galaxies with cz < 6000 kms~^, the of the TF regression drops from 3362 
to 2616 when all the modes are included, assuming a scatter of cr^ = 0.05. The rj scatter 
of the different pieces of Mark III is not identical, but the observed r] scatter after fitting 
the modes is spatially random, so a nonuniformly weighted fit will be unchanged from 
the present procedure. When a 3.5% correction is made to unbias the slope of the TF 
regression, to 7 = —0.112, the increases to 2623, which is statistically significant, but 
the resulting fiow field is negligibly changed. The derived TF slope is 2a deviant from the 
slope derived by the inverse fits for the Han-Mould sample, 7 = —0.1177 ± 0.0025 (Willick 
et al. 1995, 1996a), but this slope is derived with a very different fiow model, so we consider 
the agreement satisfactory. 

Note that the resulting per degree of freedom is 1.15, significantly greater than 
unity. The excess fiuctuations are primarily caused by small scale velocity fields and 
multi-valued zones with coherence scale below the resolution of the mode expansion. This 
is a slightly higher reduced than obtained with the mock catalog, and the increased 
chi^ could be due to slight deviations from Gaussianity in the TuUy-Fisher scatter of the 
Mark III catalog, an effect we did not include within the mock catalogs. Our pruning of 
outliers from the Mark III catalog was very modest. 

We have examined plots of the distribution of 57] in redshift and on the sky; all the 
plots are uniformly random. In Figure 7 we show the correlation function of the 77 residuals 
(see e.g. Gorski et al. 1989), before and after the modal expansion. This plot demonstrates 
the residuals to have zero coherence on scales larger than the resolution scale of the mode 
expansion. It appears as though the mode expansion has done its job; the expansion of 
Pi in a set of smooth functions has greatly reduced the scatter of the inverse TF diagram. 
The randomness of the 5r] scatter suggests that the 56 mode expansion is sufficient for 
this dataset. With a larger dataset and/or more precise distance indicators, a larger set of 
expansion coefficients would be justified. 

The resulting ITF velocity field of the Mark III galaxies is shown in slices of redshift 
space in Figure 8, and projected onto the Supergalactic plane in Figure 9. Both of these 
velocity maps are in the Local Group frame. In the nearby zone, the fiow is dominated by 
infall to Virgo (/ = 284°, 6 = 74°) and Ursa Major (145°, 65°) in the North, and modest 
outfiow from Fornax (237°, —54°) in the south. In the middle redshift zone, the dipole 
pattern persists, and is more pronounced. This is composed partially by refiex dipole of 
the motion of the LG, plus backside infall from behind the Virgo supercluster. The Hydra 
cluster seen at (270°, 26°) is observed to be falling toward us, while the foreground of 
the Centaurus region (310°, 20°) is moving away from us. In the most distant slice of 
redshift space, the dipole pattern is further enhanced, with the entire Southern galactic 
sky, including the Perseus-Pisces region (150°, —15°) and the Pavo-Indus-Telescopium 
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region (320°, —15°) flowing away from us, the Centaurus cluster flowing away rather 
substantially, but the background of Hydra and the Northern galactic cap flowing toward 
us. All of this is a combination of the reflex of the motion of the LG, plus the motion 
of the individual regions. Note the very strong shear in the observed ITF field along the 
boundary I = 300, b > 0, which divides the Hydra and Centaurus regions, and amounts to 
a change of radial peculiar velocity of approximately 900 kms""*^ in a transverse interval of 
4100 kms~^, or 20% of a Hubble flow differential motion over a substantial volume! This 
is quite a major gradient, one that appears to have been conflrmed by distance estimates 
of elliptical galaxies in Centaurus (Lynden-Bell et al. 1988, Dressier 1994). In fact the true 
gradient is limited by the angular resolution of our modal expansion, since the reported 
outflow of the Centaurus clusters is 500-2000 kms~^ in the LG frame, whereas the ITF 
measured flows for the entire Centaurus cluster are in the range 350-450 kms~^ in the LG 
frame. 

For comparison, we show in Figure 10 the IRAS predicted fleld for (3 = 0.4 on shells of 
redshift, and as projected onto the Supergalactic plane in Figure 11. Note that the Mark 
III galaxies have an IRAS predicted mean Hubble flow of 2.9%, which was removed in the 
process of smoothing the IRAS predicted field by the mode expansion. Residual fields are 
shown in Figures 12-14 for (3 = 0.4 and 0.7. Note the similarities and differences between 
these Figures and Figures 5-6 for the mock catalog. It is important to remember that 
the ITF velocity and IRAS gravity fields have been processed completely independently of 
each other. To a first approximation, the agreement between these Figures is a spectacular 
and remarkable confirmation of the consistency between the gravity and velocity fields of 
our local Universe. The IRAS galaxies must have some close relationship to the underlying 
mass distribution. The qualitative alignment of the fields is consistent with the large scale 
structure having been formed by the process of gravitational instability. 

On the other hand, the residual fields, for any value of /?, are not nearly as clean 
as seen in the mock catalogs, and are always dominated by strong dipole patterns that 
increase linearly with redshift. This is especially pronounced for the higher values of /? and 
is particularly clear in Figure 14. The IRAS maps never reproduce the strong shear seen in 
the ITF maps between the Hydra and Centaurus regions. They never reproduce the strong 
infall measured by ITF for the I = 180°, b = 45°, cz = 5000 kms~^ region. For /3 < 0.3 
(not shown), the local supercluster fiow is much too small, and the fiow toward Perseus- 
Pisces is too small. For /? = 0.4, the residual field in the local zone is still completely 
coherent, but the flow toward Perseus-Pisces and the average of Hydra and Centaurus is 
approximately correct. For (3 = 0.7 the relative motion toward Virgo matches the ITF 
measurement (—340 kms~^ in the LG frame) but the flow away from Perseus-Pisces and 
toward the Great Attractor complex is now too large. 

4-3 The Dipole Residual on Shells 

The residual fleld between the ITF and IRAS flelds is clearly dominated by the dipole 
component. ND94 emphasize that the dipole component of the radial peculiar velocity fleld 
on a shell, relative to the LG, is generated entirely by the mass distribution within that 
shell; the exterior mass distribution affects only quadrupole and higher I perturbations. 
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Thus it is critical to seek agreement of the dipole components of the fields on modest scales, 
cz < 6000 kms~^, where the data are best. 

Although the large spatial inhomogeneity of the Mark III data strongly mixes the 
modes together, it is possible to extract the dipole component of the field. Consider 
binning the Mark III data into shells of redshift 1000 kms~^ thick. Then compute the 
three components of a dipole fD by summing the radial peculiar velocity Vj over the j 
galaxies in the shell: 

EjVjWjCOs{e) 

where $ is the angle from the three unit vectors x, y, z- By using a weight function Wj 
given by (5^4)"^, where 864 is the angle of galaxy j to its fourth nearest neighbor on the 
given shell, the computation of vd is approximately equal area weighted, as is appropriate 
for comparison to the original functional expansion. 

An example is shown in Figure 15, where we plot the ITF and IRAS (for f3 = 0.6) 
derived dipoles in six shells within 6000 kms~^. The open and closed symbols are derived 
from the irregularly spaced Mark III galaxies, while the solid lines are original IRAS dipole 
components, for /? = 0.6, based on the full 1.2 Jy sample. Note that the IRAS derived 
points, drawn from a less than optimal sky distribution, still match the smooth curve 
very well, demonstrating that recovery of the dipole from the irregularly spaced Mark III 
sample is not a problem. The open squares showing the ITF dipole are consistent with 
the IRAS dipole in the Z component, and the total amplitude is a good match as well. 
But the X and Y components are seriously discrepant. The X component is consistent for 
the first three shells, but then deviates substantially, while the Y component is discrepant 
for all shells; IRAS predicts a substantial reflex dipole in the Y component that should 
grow with redshift, but little is detected in the ITF data. For those who have followed 
this subject for some time, this discrepancy is not new (see e.g. the comparison of Faber 
and Burstein (1988) with the discussion by Strauss and Davis (1988) and Yahil (1988)). 
The IRAS field expects to find substantial "backside infall" at the Centaurus region and 
beyond, but such a signature is very weak in the Mark III data. We emphasize once again 
that the comparison under consideration here is for spiral galaxies only; the discrepancy 
is observed for both spiral and elliptical galaxies. 

5. Likelihood Analysis of Modal Comparison 

5. 1 The scatter of ITF and IRAS coefficients 

Given the measurement of the coefficients derived from the Mark III and IRAS fields, 
how can we best compare them? In the discussion above, we compared the fields galaxy by 
galaxy, and by means of the area weighted dipole on shells. The sky maps are appropriate 
for visualization; however the fields as seen in space are highly correlated because of the 
strong covariance of the errors in the ITF velocities as well as in the IRAS gravity. Since the 
number of modes which describe the fields is significantly less than the number of galaxies 
in the TuUy-Fisher sample, a galaxy by galaxy comparison requires inversion of an unduly 
large covariance matrix. Therefore we work directly with the recovered mode coefficients 
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from the ITF and the IRAS gravity in order to assess the quantitative agreement between 
the velocity and gravity fields. 

Figure 16 plots the coefficients atf versus CKj^os for the real data, using /? = 0.6 for the 
IRAS coefficients. Also shown are the coefficients for one realization of the mock catalog; 
in this case we can also plot the coefficients atf versus atme- Note in the mock case that 
the coefficients track the diagonal line extremely well, consistent with (3 — 1. As expected, 
the noise is increased for the ITF versus IRAS comparison. Different realizations of the 
mock catalog all have consistent behavior in their mode-mode scatter diagrams. 

5.2 The error matrix on ITF and IRAS coefficients 

As is shown in ND, the orthonormality condition leads to the following simple form 
for the error matrix of the ITF coefficients, namely. 



{6al, 6a{,) = 6^., . (20) 
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In Figure (16), the plot of at/ versus atme for the mock catalog has noise that is completely 
consistent with equation (20), as expected. 

For the Mark III catalog, six separate samples, with differing photometric bandpasses 
and velocity width methods, have been pieced together in the manner described above. 
The offsets in the magnitudes and linewidths for this adjustment have uncertainty in the 
mean, and this will translate into a non-diagonal contribution to the covariance matrix 
< 5a-[j Salj- >. We estimated this contribution by generating 20 realizations of the Mark 
III catalog with gaussian errors in the zero point of the different pieces appropriate for each 
subsection. We then directly estimated the covariance matrix generated by these offsets, 
and included them in the analysis below. In fact, the errors in the zero points are small 
and make no difference to the error budget. 

The situation for the modes of the IRAS predicted velocity field is not so simple, 
however. The error matrix < Sal^^^ ^'^Iras contructed from the covariance 

matrix of the IRAS predicted velocities of galaxies. Using Equation (15), and linearizing 
the errors in the velocities about the predicted values yields, 

-< ««^r-a. ««^r-a. > ^ ) (1 - (z,,/^,)) (1 - («,,/^,0) ' * 

(21) 

where < 5uiSui' > is the covariance matrix of the IRAS predicted velocities of galaxies j 
and j' after subtracting the Hubble-like flow as described in 2.2.2. Note that even if this 
velocity covariance matrix is diagonal with equal elements, the covariance matrix of the 
coefficients will not be diagonal, as the modes F are orthogonal with unit weighting of 
each point, whereas the sum above is weighted approximately as z~'^ . 

The IRAS velocity field is a very nonlocal quantity, effectively computed by counting 
the vector sum of all IRAS galaxies around a given point, weighted by the inverse of the 
square of the real space separation and by the inverse of the selection function. This 
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sum is computed over a discrete set of points in a finite catalog, and therefore must 
suffer from Poisson noise and from variance induced by the assumption that the mass 
per IRAS galaxy is assumed constant for each redshift shell (c.f. Strauss et al. 1992). In 
addition, the computed velocities assume linear theory, yet there exist high density zones 
within the galaxy distribution for which linear theory is known to be inaccurate and for 
which the deviations are likely to be coherent (as in triple-valued zones or zones which 
have shell crossed). There are incoherent small-scale nonlinear velocities unresolved by the 
IRAS procedure which add noise as well. Furthermore, the IRAS predicted velocities are 
sensitive to the assumed transformation from heliocentric to LG frame velocities, which 
might be in error by perhaps 50-100 kms~^. All of these error sources contribute to errors 
in the IRAS reconstructed velocities, and most of these errors are coherent. 

We can express the IRAS velocity covariance as a sum over these terms, writing 

< SuiSuj >= CsN{i,j) + CNL{i,j) + f^iGCOs(%) (22) 

where CsN{hj) is the covariance in the predicted peculiar velocities due to the Poisson 
noise, Cnl is a term to describe the deviations from linear theory, and is the uncer- 
tainty in the transformation to LG frame. We shall set (Jlg = 100 kms~^. Note that this 
term alone generates complete covariance of all IRAS predicted peculiar velocities. 

The Poisson noise term CsNihj) can easily be derived from a distribution of points 
in real space (Yahil et al. 1991) where, according to linear theory, the peculiar velocity 
field is related to the distribution of galaxies by the Poisson equation, and the variance 
CsN would scale as /?^. However, the galaxies are observed in redshift space and the error 
in the predicted velocity field must be computed according to the relation (4), for which 
the scaling with /? is not trivial. An analytical expression for Csn for the velocities from 
(4) is cumbersome and numerically expensive to evaluate. Thus we choose to compute 
the covariance matrix Csn by generating 18 bootstrap realizations of the observed 1.2 
Jy IRAS galaxy distribution by replacing each observed galaxy with a number of points 
drawn from a Poisson distribution with mean of unity. Then we compute the velocity fields 
from these realizations by the same algorithm used in the derivation of the velocity field 
from the observed galaxy distribution. For galaxy i of the Mark III sample, we tabulate 
the differences, 5ui, between the velocity obtained from each of 18 bootstrap realizations 
and the velocity as predicted from the actual distribution of galaxies, and evaluate the 
covariance matrix Csn by averaging the product SuiSuj from all the bootstrap realizations. 
This process is computed for several different values of /3. 

The term Cjvl(^, j) describing nonlinear effects can be modelled as 

CNLiiJ) = olb% + Hsy^ii) + Sy,{j)f CXp [- > (23) 

where particles have an assumed incoherent small scale 1-D velocity dispersion a^-, but in 
regions of high shear can also display coherent errors, with coherence length Ucoh which we 
set equal to 300 kms~^. Miller et al. (1996) have directly measured (7„ for the IRAS 1.2 
Jy survey, and derive a best value cr^ = 137 kms~^. Willick et al. (1996c) find a best value 
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of (T„ ~ 150 kms"-*^ for the Mark III galaxies using the VELMOD analysis (i.e. "Method 
11"*"" as described by Strauss and Willick (1995)). The Mark III and IRAS galaxies need 
not have the same thermal dispersion, but the agreement is comforting. We shall assume 
a value 0"^ = 150 kms~^. The coherent error we assume is proportional to the square of 
the average shear in the IRAS derived velocity field, Shii) = 100 min(l, \dvp/dz\) kms"-*^, 
which is an approximately linear function of /3. This term is designed to deweight points 
in regions of probable shell crossing, where the IRAS predicted shear is large. 

5.3 Pseudo analysis 

Our basic assumption is that, for the true (3, the ITF and the IRAS velocity fields are 
noisy realizations of the true underlying velocity field. Let be the coefficients obtained 
from the modal expansion of the true underlying field. Each estimated coefficient is derived 
from a sum over a large set of numbers, so by the central limit theorem we expect each 
estimated coefficent to have a normal distribution, even if the underlying velocity field 
or TuUy-Fisher data has nongaussian scatter. Since the Mark III and the IRAS data are 
independent we write the probability for measuring a^^^^ and a^j as. 



Aexp 



(X?/ + Xlas) 



(24) 



where we have defined. 



'tf 



and. 



Xiras 



= - «iaJ(M,y )(«^' - (4ras)- 



3,3 



Note that x^f does not depend on Given /3, we seek the true coefficients, a-' , and /3, 
by maximizing the probability (24). Therefore our best estimates are obtained by solving. 



dP/da^ = 0, 



(25) 



and. 



dP/dl3 = 0. 



(26) 



Since the normalization of the probability (24) does not depend on the coefficients a-' , the 
condition (25) implies that the parameters which maximize P for a given j3 renders a 
mimimum in the total = xHras + X?/- However, since the error matrix M and hence 
the normalization factor A strongly depend on /3, the maximum of P with respect to (3 
docs not necessary coincide with the minimum of x^. Though, at a given /3, the condition 
(25) yields a set of linear equations for which can be solved using standard methods, 
the condition (26) leads, in general, to a complicated nonlinear algebraic equation for j3. 
Therefore, a simultaneous estimation of the coefficients and /5 is a problem of nonlinear 
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minimization. Fortunately, in this particular case, it is easy to see that solving (25) and 

(26) simultaneously can be achieved by the following scheme: solve (25) for at a given 
f3 and compute the corresponding value of P{P)- The estimate for the true /3 is the value 
which renders P{(3) a maximum. The condition (25) yields, 

a ■ (la-/ + M"^) =ottf la-/ + a^ras ■ M'K (27) 

By substituting (27) for in the total x^, we have, 

X^\dxydaJ=0 Xtf + Xiras = J2 ~ ^iras) (1^1/ + M)~! (a^j - a^^^) , (28) 

3,j' 

which reduces to the pseudo x^ suggested by Press et al. (1992) for problems with errors 
in both the x and y directions. Note that the difference between the values of P which 
minimize the x^(/3) and maximize the probability P{P) must be consistent with the errors 
in their estimates provided that the ITF and IRAS coefficients differ only because of the 
noise in their measurements. 

Comparing the TuUy-Fisher coefficients to IRAS derived coefficients a^^as for the 
mock catalog gives — 53.5 for 56 degrees of freedom. Thus the fit of the coefficients 
of the mock catalogs is entirely consistent with the expected noise, and this consistency is 
seen in all the realizations of the mock catalog. When we apply equation (22) to the mock 
catalog, comparing the TuUy Fisher coefficients to the true coefficients alme^ while 
holding Mij nonzero, we obtain = 30.2 for 56 degrees of freedom. This low results 
from inappropriate inclusion of the IRAS like covariance; the approximate halving of the 
demonstrates that the errors introduced by the imprecision of the IRAS estimates are 
similar in magnitude to those of the ITF estimate. 

The results of applying this same procedure to the real data is given in Table 1, where 
we list x^ versus For each value of /? we first computed M^^ (equation 16) and then 
evaluated equation (23). The table shows that the minimum = 100 (55 d.o.f.) at 
/3 = 0.5. Given that x^ for u degrees of freedom has an expected mean of u and a variance 
of 2iy, the best x^ is 4.3a" deviant from its expected value. The fit is statistically very 
improbable and can be rejected. That is, although the mock catalogs demonstrate that 
the procedure is capable of producing very acceptable x^ results, the actual comparison of 
Mark III to IRAS leads to the conclusion that the fields are statistically inconsistent with 
each other. The details of the best /3 and the level of inconsistency are weak functions 
of how the catalogs are limited; we would derive a higher (3 if we limited the catalogs to 
cz < 4000 kms~-^, but the fit would remain unacceptable. 

Because the residual field is dominated by a strong dipole and because we believe we 
understand the error matrix, it is not permissible to simply double the covariance matrix 
to bring the x^ values to reasonable levels. Thus we urge extreme caution in concluding 
that /3 = 0.5 ± 0.1 is in any sense a "best fit" solution. 

6. Discussion 
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6.1 Alternative Methods and Results 

The ITF analysis is only one of several, independent methods of comparing the TuUy- 
Fisher catalogs with the redshift survey catalogs. Given the difficulty of the comparison 

and the subtlety of the error analysis, it is worthwhile to investigate alternative approaches 
and to expect consistent results. Only then will we have confidence in the derived value of 
(3. 

The great advantage of the ITF scheme is that it allows the comparison of the velocity 
and gravity field mode by mode, with identical resolution to each field. The Monte-Carlo 
simulations with realistic mock catalogs confirm that the bias in the extracted velocity field 
is negligible. The method allows for the gradient in radial resolution of the two fields, so it 
is capable of extracting the maximum useful information for each field. One can rank order 
the modes by their significance, or contribution to x^, keeping only the most significant 
modes, so as to provide a very compact description of the field. This is an example of a 
Karhunen-Loeve transformation (Therrien 1992). Recent applications of this technique to 
astrophysics include the CMBR fluctuations of the COBE data (Gorski 1994, Bunn 1995, 
Bond 1995), the power spectral analysis of fluctuations in redshift surveys (Vogeley and 
Szalay 1996), and the spectral classiflcation of galaxies (Connolly et al. 1995). 

The "VELMOD" approach of Willick et al. (1996c; cf. the preliminary discussion 
presented by Strauss & Willick 1995, § 8.1.3) is in several respects complementary to the 
analysis presented here. VELMOD also makes a comparison of the Mark III TF data 
and the gravity field derived from the IRAS 1.2Jy redshift survey. It differs from the 
ITF analysis in explicitly accounting for the small-scale velocity dispersion o"„, which it 
treats as a free parameter, and in allowing for the possibility of multivalued redshift- 
distance diagrams; the ITF reconstruction implicitly assumes a unique redshift-distance 
mapping. Implementation of VELMOD is done using grid-based, real-space IRAS velocity 
reconstructions, as opposed to the noniterative redshift-space IRAS velocity fields used 
here. Thus, VELMOD and ITF do not use identical gravity fields, although tests have 
shown them to be very similar. 

There are several advantages to VELMOD . First, in accounting for small-scale veloc- 
ities and arbitrary redshift-distance diagrams, it determines the probabilities of individual 
data points without some of the approximations used here, and thus provides a reliable 
maximum likelihood statistic. It does not require that the input TF catalog be placed on 
a common system, but instead treats the individual sample TF relations as free parame- 
ters. This makes it free of the uncertainties induced by the overlap-comparison procedure 
described above (§ 4.1). Finally, VELMOD is particularly well-suited to analysis of the 
nonlinear regime, which holds the promise of breaking the fl-b degeneracy. However, 
VELMOD has an important weakness from which neither ITF nor POTENT suffers: it 
does not permit the construction of a velocity field from the TF data alone; instead, it 
assumes that the IRAS velocity field, for some value of (3 (and 6), is a good fit to the data, 
and selects the best (3 by maximizing likelihood. Yet as we have shown in this paper, the 
IRAS velocity field does not fit the data well throughout the sampled volume. In such 
a circumstance, the meaning of the maximum likelihood statistic is unclear. In addition, 
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VELMOD is not well suited to a proper analysis of covariant velocity field prediction errors, 
such as we have carried out here. 

A final distinction derives from the depth of the respective analyses. VELMOD is 

rather cpu-intensive and, moreover, differs meaningfully from analyses such as ITF only 
when au is a substantial fraction of the TF distance error. The analysis of Willick et al. 
(1996c) is thus restricted to redshifts cz < 3000 kms~^. In contrast, we have extended 
our computationally more efficient analysis to twice that distance. As noted above, our 
low indicated value of /3 = 0.4 — 0.6 is in part a consequence of discrepancies between the 
ITF velocity field and the IRAS gravity field that are manifest principally at distances 
^ 4000 kms~^. Willick et al. (1996c) analyze their likelihood fit residuals and find them 
to be nearly incoherent across the sky, which refiects the fairly good agreement between 
IRAS and the TF data within ~ 3000 kms"-^. Preliminary results seem consistent with 
our "best" values of 

Both the POTENT and ITF methods have been validated by extensive testing with 
mock catalogs, so it seems reasonable that they should yield consistent results. We here 
list several several relative advantages and disadvantages between the ITF method and the 
POTENT/IRAS comparison. 

(1) POTENT provides a measure of the divergence of the velocity field which can be 
compared to the directly observed galaxy density field. Both fields are local, an 
obvious advantage. The ITF field must be compared to the non-local IRAS derived 
velocity field. These fields are both susceptible to systematic errors within the data. 

(2) POTENT required a large, uniform smoothing to properly define the smoothed veloc- 
ity and velocity potential. Dekel et al. (1993) have been using a gaussian smoothing 
of width a — 1200 kms~^, sufficient to filter out the infall of the Local Group toward 
the Virgo cluster. The ITF method can work with much smaller smoothing, and the 
smoothing length can be dependent on redshift. This is quite advantageous, since 
the statistical precision of both the velocity and gravity field are decreasing functions 
of distance. On the other hand, the ITF method can readily accomodate uniform 
smoothing, if that is desired. 

(3) The POTENT analysis is dependent on a calibrated forward TuUy-Fisher relation (i.e. 
M = M{ri), where M is the absolute magnitude and ry = log(At;) is the velocity width 
parameter), and the analysis is sensitive to Malmquist bias (see Strauss and Willick 
1995 for details). The ITF analysis uses an inverse TF relation (i.e. r] — r]{M)); the 
parameters of this relation can be fitted at the same time as the mode coefficients. As 
long as the galaxy selection is independent of linewidth, this method does not suffer 
from Malmquist bias. The ITF method is biased in the presence of a thermally hot 
velocity field, where multistreaming is very pronoimced; however, given the estimated 
amplitude of the small scale peculiar velocity field of galaxies, this bias is quite modest, 
as described above (see ND95 for details). 

(4) The ITF method is simply a redshift space smoothing algorithm for the TuUy-Fisher 
data that enables one to get a measure of the velocity field without binning the 
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data. It is not intended for the construction of full-sky velocity maps, and such maps 
are not needed for comparison to the gravity field. The method compares the two 
fields only where there are data. By avoiding binning of the data, it eliminates the 
"sampling gradient bias" that is a serious problem for POTENT. Both the ITF and 
IRAS analyses are here carried out in observed redshift space, whereas most previous 
analyses, including POTENT, have been done in real space. Since the inversion to 
real space is not unique in the multi-valued zones, attempting to perform the ITF 
analysis in real space would present additional very serious complications. However, 
with sufficiently precise distance indicators, working in redshift space will not take 
full advantage of the data, since multivalued zones could be resolved (e.g. Strauss and 
Willick 1995). 

(5) By filtering the density and velocity fields with identical functions, one can be con- 
fident the fields have identical resolution. In the POTENT analysis, assessing the 
relative resolution of the galaxy density field compared to the inferred mass density 
field is not straightforward. Previous velocity-velocity comparisons have been espe- 
cially compromised by the inability to filter the velocity field. Although we have here 
chosen a very general set of functions for the mode expansion, if desired, it would be 
straightforward to supplement this expansion with additional functions to describe 
known features, such as Virgocentric infall. This would be an appropriate route for 
datasets smaller or more concentrated in selected directions than Mark III. 

(6) The full covariance matrix of the coefficients of the ITF and IRAS fields can be com- 
puted, and because the matrix is modest in size, it can be readily inverted. This 
is essential for a proper statistical assessment of the fit. POTENT compares veloc- 
ity/dentiy fields computed on a grid, and the error in P must be estimated using 
sophisticated Monte Carlos methods. All previous point by point velocity-velocity 
comparisons have failed to include this essential covariance in the error matrix of the 
gravity field. 



6.2: Why Don't the Fields Agree? 

The ITF method has led to an unambiguous measure of the TuUy-Fisher velocity field, 
one we feel is correct within the limited resolution of our chosen expansion. The biases 
generated by small-scale thermal velocities are modest. There is undoubtedly additional 
information that could be gleaned from the Mark III database, as our basis function of 
modes is not necessarily optimal for the derivation of the divergence of the flows, the 
critical information extracted in the POTENT analysis. The tests with the mock catalog 
validate our approach, and demonstrate that comparison of the Mark III data with the 
IRAS derived field should yield consistent fields. The error matrix of both the Mark III 
and IRAS derived mode coefficients is, we believe, well understood. Yet the comparison 
of the modes conclusively demonstrates that they are not consistent within their errors 
at the 4-0" level. We have already used a very generous TuUy-Fisher error, cr^ = 0.05 
(corresponding to ctm ~ 0.45). Increasing this error will reduce the formal but will not 
remove the coherent residuals. 



25 



The nature of the inconsistency is dependent on the depth of the sample under con- 
sideration. To fit the Local Supercluster will require /3 ^ 0.7, but the large dipole residual 
at large scales demands a smaller /3. Thus it is not surprising that the "best" /3 for our 
analysis, (3 — 0.4 — 6, is slightly lower than suggested by VELMOD , (3 — 0.57, as that 
analysis is confined to galaxies with cz < 3000 kms~^. 

The 4-(j discrepancy between the IRAS derived velocity field and the Mark III velocity 
field suggests that the problem is not due to random errors in either field, but to systematic 
problems in one or both datasets. There are many possible sources of systematic errors, 
and we here list and discuss the most obvious among them. It is very likely that some 
combination of these errors is responsible for the statistically poor correspondence of the 
fields, and until we can deal with the major problems, all derivations of (3 from fiow analysis 
must be considered suspect. 

(1.) The IRAS field is missing the early type galaxies, so it ignores cluster centers. Careful 
comparison of the counts of IRAS galaxies to optically selected galaxies indeed shows 
an underrepresenation of the galaxies within the centers of clusters (e.g. Virgo's core 
is underrepresented by a factor of 2). However the cores of clusters do not dominate 
the mass of the superclusters (recall that the cluster mass typically diverges as radius 
R), and the IRAS galaxies appear to be good tracers of the general field population 
(Strauss et al. 1992). "Boosting" known cluster centers to compensate for the missing 
early type galaxies does not lead to a better fit between the gravity and velocity fields, 
but it will reduce the best estimate of /? by a factor of ft! 1.3. 

(2.) The IRAS field does not have whole sky coverage and is missing a prominent cluster 
in the Zone of Avoidance. Recently Kraan-Korteweg et al. (1996) have reported on a 
Coma-like cluster, A3627, with estimated mass 2.5 x lO^^/i-i Mq at position 1=325°, 
b=— 7°, czlg ~ 4700 kms"-*^. Although the IRAS survey does include galaxies to 
\b\ > 5°, and the observed galaxy density is interpolated through the Galactic plane, 
a cluster of entirely early type galaxies would have been severely undercounted. The 
IRAS catalog has an overdensity of galaxies at the position of A3627, but it is not 
sufficient to comprise a large Coma-like cluster, and it is not by any means the most 
massive IRAS cluster in this vicinity. Could this possibly explain the strong shear ob- 
served in the ITF analysis between the Hydra and Centaurus regions? Suppose that 
the IRAS survey has given no mass to A3627. Then for the parameters above, the ne- 
glect of A3627 has led to an error of the peculiar motion of the LG of ~ SOJl" "^ km s~^. 
We can readily compute the missing differential motion between the Centaurus clus- 
ters at (l,b) = (302°, 22°) and the LG caused by A3627. If Centaurus is at distance 
of ft 3800 kms~^, as suggested by the IRAS field, then its acceleration toward A3627 
is in the transverse direction. If the distance of the Centaurus clusters is lowered 
to 2500 kms~^, as suggested by several distance estimates (Lyndcn-Bell et al. 1988; 
Dressier 1994), then the differential radial acceleration between Centaurus and the LG 
induced by A3627 is only 4Q,~-^ kms~-^. The differential influence of A3627 on Hydra 
(at (1, b, cz) = (275°, 20°, 4500 km/s)) and the LG is approximately 75fl-°-^ kms-^ 
These results scale linearly with the mass of A3627. Thus as one increases the mass 
of A3627, the residuals toward Hydra become more negative whereas the residuals 
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toward Centaurus are changed by a negligible amount. This is the wrong sense, since 
the residuals toward Centaurus are far worse than those toward Hydra. Furthermore, 
the effect of a missing cluster like A3627 is modest compared to the residuals to be 
explained. Examination of the IRAS density field at this redshift demonstrates that 
the total mass overdensity in the Hydra-Centaurus-Pavo-Indus complex greatly ex- 
ceeds that of any missing cluster center, and that the addition of the core of A3627 
makes no difference to the overall density field. 

(3.) The IRAS field has poor convergence to the CMB dipole either because of dilute sam- 
pling noise at large distance or because of an improper transformation to the Local 
Group frame. The ITF method is sensitive to the choice of the LG frame. The Mark 
III data, when analysed by the POTENT algorithm, exhibits a bulk motion of 250- 
300 kms~^ that is not observed in the IRAS derived gravity field (Dekel 1995). If 
this is due to an incorrect choice of the LG frame, then the effect is in the direction of 
the nearby Maffei/Dwingeloo group of galaxies (Kraan-Korteweg et al. 1994), which 
are highly obscured and are not included in the IRAS gravity field. If they are the 
source of this anomalous motion, then they should be blueshifted by our infall toward 
them. However they are observed to have redshifts, not blueshifts, and so must move 
coherently with the Milky Way. Beyond the local galaxies in this direction is a void, 
which should repel us, not attract us. It is difficult to understand how the LG could 
have a local source of acceleration leading to a velocity as large as that induced by 
the Virgo supercluster, without leaving any signature in the local galaxy distribution. 
The local flow field has only modest shear, and a large error in the LG frame definition 
is unlikely. This is apparent in Figures 8, 14, and 15; for shells within 2000 kms~^, 
the dipole residual between the ITF and IRAS fields is very modest. An error in the 
LG frame would lead to dipole residuals constant with redshift, whereas the observed 
dipole residual grows approximately linearly with distance (Figures 14 and 15). The 
local flow is rather quiet (Peebles 1992; Schlegel et al. 1994) and shows no evidence 
of the very large local shear that would result if the LG frame were offset by several 
hundred kms~^. 

As a test of the redefinition of the LG, we have generated IRAS fields and ITF fields 
with the LG frame offset toward Maffei by 200 kms~^. Such an offset leads to smaller 
residuals toward Perseus-Pisces for high values of /3, compared to Figure 14, but 
instead there are coherent errors toward Hydra and other nearby regions, as expected. 
We were unable to achieve a better fit by adjusting the LG frame. 

The diluteness of the IRAS density field could lead to substantial errors in the pre- 
dicted Uiras field, especially from matter at redshifts cz > 6000 kms~^ (Davis et al. 
1991). The mock catalogs demonstrate the degree of expected error. However, errors 
in the density field of external matter will not affect the / = 1 residuals in the LG 
frame. These dipole residuals, if real, are indicative of a mass field not traced by the 
IRAS galaxies. Since we expect linear theory to apply, we can simply estimate the 
unseen mass density by taking the divergence of the residual dipole field. Making 
allowance for the transverse fields by assuming potential fiow, and noting that the 
residual field scales approximately linearly with redshift, we find a residual mass field 
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5p oc cos^, where d is the angle from the apex direction. This is a very strange density 
field, one that is constant in redshift but tilted toward a particular direction. Thus the 
residual dipole field cannot be considered physically reasonable; it must be indicative 
of some systematic discrepancy. 

The reason our results are not consistent with the high (3 values derived from the 
density /density comparisons (Dekel et al. 1993) using the POTENT algorithm is clear 
from Figure 14. For large values of j3 the dipole residual growing linearly with z is very 
pronounced. It cannot be removed by a transformation of the LG center, nor can it 
be explained by an external mass field. It can only be reduced by elimination of some 
unknown systematic error in one of the catalogs, or the reduction of (3. The POTENT 
algorithm, by taking the derivatives of the Mark III velocity field, is less sensitive to 
this large scale dipole, but it is more sensitive to higher spatial frequency components 
of the flow fleld. This very different weighting of the data, combined with the poor 
quality of the flt of IRAS to Mark III for any value of /3, is sufficient to explain why 
it is not unreasonable for the ITF/VELMOD algorithms to yield "best" values that 
are inconsistent with the best estimates derived from POTENT using the same data. 
One could argue that the difference in the derived values of (3 is the result of a strong 
scale dependence of the biasing relation. However, as shown by Kauffman, Nusser and 
Steinmetz (1996), all physically reasonable bias mechanisms for hierarchical models 
of structure formation are well approximated as linear, scale independent biasing for 
scales larger than a few Mpc. 

) The Mark III database is not truly uniform; the effort to unify the five individual 
sections of Mark III to a common linewidth, magnitude system is somehow deficient. 
The Mark III sample is a collection of available data and samples different directions 
far from uniformly. There is a large concentration of galaxies toward Perseus-Pisces 
contributed by Willick (1990), and a very large concentration of galaxies in the South- 
ern sky contributed by Mathewson et al. (1992) on the opposite side of the sky. More 
importantly, different techniques are used in differet directions. Any inconsistencies 
within the Mark III catalog are likely to translate into dipole effects, and indeed the 
dominant residuals between Mark III and IRAS is a dipole aligned along the center 
of the MAT subsample. 

Willick et al. (1996a,b) describe the comparison of the different pieces of the Mark III 
sample. It is obviously critically important to insure that the relative photometry and 
linewidth measurements are standardized, and it is possible that the Mark III data has 
zero point calibration problems between its different subpieces. As described above, we 
have tested the sensitivity of the ITF field to the known statistical uncertainty of the 
zero-point offsets, finding a negligible effect. If calibration mismatches are the source of 
the problem, then the error must be considerably in excess of the measured statistical 
errors. It would be especially valuable to obtain additional data to better tie together 
the different pieces of Mark III, and to improve the sky coverage in the North and away 
from the dense cluster centers. The Mark III sample is vulnerable to any systematic 
nonuniformities between Han-Mould Northern versus Southern clusters, because this 
subsample is used to renormalize the different pieces of Mark III to one system. A new 



28 



survey of TuUy- Fisher distances for 300 galaxies in the zone 5000 < cz < 7000 kms~^ 
around the entire sky is just now getting under way (Strauss et al. 1996), with the 
goal of better calibration of the whole-sky samples. 

A considerable portion of the Mark III sample is derived from relatively low galac- 
tic latitude data, where corrections for galactic extinction are significant. Recently 
Schlegel, Finkbeiner, and Davis (1996) have been studying whether a new extinction 
map, based on the COBE/DIRBE data, leads to significantly different flow fields. The 
question is still under investigation, but the answer appears to be no. 

(5.) The Tully-Fisher calibration is dependent on environment and the residual field is 
caused by changes in the TF parameters for different sections of the sky. Although it is 
impossible to prove that the TF relation is independent of environment, it is not easy to 
imagine schemes whereby the relation is coherently different in one portion of the sky 
compared to another. It is much easier to imagine environmental dependence leading 
to larger scatter of the TF relation. It is important to note that the most serious zone 
of discrepancy, in the Centaurus region, has been extensively studied with the — a 
relation for ellipticals (Lynden-Bell et al. 1988) as well as recent surface brightness 
fluctuation measurements using ellipticals (Dressier 1994; Tonry and Dressier 1995). 
All three distance indicators suggest strong outflow in the Centaurus region, which 
the IRAS maps fail to predict. 

Improving the IRAS maps by going to denser sampling, either using the PSC catalog 
(Saunders et al. 1995), or by supplementing the IRAS catalog with optical data in regions 
not too badly extincted (e.g. Santiago et al. 1995), will be an important next step to 
attempt to resolve the discrepancies here presented. Based on previous analyses (Hudson 
1994; Freudling & da Costa 1994), these improvements in the gravity map are unlikely to 
resolve the problem. More uniform full-sky Tully-Fisher catalogs, where the line widths 
and magnitudes are estimated in a consistent fashion for the entire sample, will be an 
enormous step forward, but this is a major undertaking. 

7. Summary 

We have presented a method for comparing the velocity fleld derived from Tully-Fisher 
type data with the gravity fleld derived from full-sky redshift surveys of galaxies, such as 
the IRAS surveys. The chief strength of the method is that the two flelds are flltered 
through the same set of low resolution modes, so that one is assured they indeed have 
the same resolution. In spite of working in the 'inverse' direction, we are able to produce 
pictures of the measured flelds. 

Although the method works extremely well in mock catalogs, with modal coefficients 
fully consistent with the expected noise, we flnd a much poorer agreement for the real 
data. This discrepancy is not consistent with random errors, and may be indicative of 
a systematic error in one or both of the datasets. We derive a "best" value of /3 = 
0.4 — 0.6, in reasonable agreement with the value derived from the VELMOD analysis, 
though considerably smaller than the values (/3 ~ 1.0) typically obtained from applications 
of the POTENT algorithm. Because the for the flt of the ITF velocity fleld to the 
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IRAS gravity field is 100 for 55 degrees of freedom, the fields cannot be said to agree 
in detail. We therefore urge extreme caution in concluding that /3 has been measured in 
large scale flows. It is worth emphasizing that there is qualitative agreement between the 
velocity and gravity fields, particularly at distances ^ 3000 kms~^. However, until the 
discrepancies on larger scales can be resolved, and all methods give consistent answers, we 
must prudently consider the value of /? to be an open question. 

Future work should greatly improve this situation, as more data, and more precise 
data, become available. The Mark III catalog is not optimal for the analysis presented 
here, since the sky coverage is relatively nonuniform. For a given allocation of telescope 
time, the optimal strategy is to sample the sky to a given depth as uniformly as possible. 
Full sky samples are difficult to calibrate in a uniform manner, but are worth the effort, 
since they provide the only route to estimation of the low order multipoles of the local 
velocity field. 
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Figure 1: The amplitude of tlie radial variation of the mode velocities versus cz 
for 1=0-3. Note that the radial resolution drops sharply with redshift. 



35 



this figure available by anonymous ftp at magicbean.berkeley.edu/pub/marc/m3itf/ 



Figure 2: The Uuf sky projection as seen in the LG frame for the mock catalog, 
in galactic coordinates. The open circles are points that are flowing inward and 
the stars are points flowing outward. The size of the symbols is proportional to 
the flow velocity, with a key showing 500 kms""*^ flow. Note the three separate 
panels of redshift interval. Note also the dominance of the dipole mode, especially 
for the outer shell. This is the signature of the reflex motion of the LG. 
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Figure 3: The unf field for points witliin 30° of the Supergalactic Plane, plotted 
in this plane, for the mock catalog, in the LG frame. Solid symbols and solid lines 
are outflowing "galaxies" with the symbol drawn at the distance and the line 
drawn to the measured LG redshift. Open symbols and dashed lines correspond 
to "galaxies" infalling to the LG. Note again the strong dipole pattern of the flow 
at ^ > 3000 kms-^ 
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Figure 4: (a) The true, unfiltered peculiar velocities, Utme of the points in the 

mock catalog, plotted in the Supergalactic Plane, (b:) The true peculiar velocities, 
J- 

Utrue of the mock catalog points, after filtering by the 56 mode expansion. Note 
how the large scale behavior of panel a is reproduced in panel b. The high spatial 
frequency signal, which is very difficult to reproduce in the IRAS gravity field, is 
effectively eliminated in Figure 4b. (c:) The IRAS predicted velocity field for the 
mock catalog, Uiras, plotted in the Supgggalactic Plane, (d:) The same as panel 
(c), but now filtered by the 56 mode expansion, uj^^^. 




Figure 5: (a:) The residuals uuf — ui^ue ^ function of cz of points in the mock 
catalog within 60° of the apex and anti-apex direction (l,b) = (270,30), before 
correction for the bias in the slope parameter 7. (b:) The same plot unf — w/^^e 

after the bias in the slope parameter is removed, (c:) The residuals of unf —uf^as 
for the same galaxies. 
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Figure 6a: Residuals uuf — ui^^^^, in sky slices for mock catalog, using the same 
notation as in Figure 2. 
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Figure 6b: Residuals uuf — u^^^g in the same projection. Note the small 
amplitude of these residuals versus Figure 2, and the absence of significant dipole 
contribution to these residuals. Note that the amplitude of the residuals in (6b) 
are somewhat larger than those of (6a), because of the IRAS contributed noise. 
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Figure 7: The auto-correlation function of the ry residuals of the Mark III galaxies 
versus redshift space separation. The dashed curve results when the 56 mode 
coefficients are set to 0, while the solid curve is after the 56 modes are fit to the 
Mark III data. 
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Figure 8: The ITF velocity field u^f for the Mark III galaxies using the same 
sky projection as Figure 2. 
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Figure 9: The ITF velocity field Unf: for Mark III galaxies within 30° of the 
supergalactic plane, using the same notation as in Figure 3. 
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Figure 10: The sky projection of the IRAS predicted flow fleld w^ras for /3 = 0.4 
for the Mark III galaxies. 
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Figure 11: The supergalactic plane projection for the /i2v4iS' predicted flow field 
iJ'iras for P = 0-4 for the Mark III galaxies. 
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Figure 12: The sky projection of the residuals Uuf — Uiras for (3 = 0.4. This is 
simply the difference between Figures 8 and 10. Compare this field to Figure 6b. 
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Figure 13: The same as Figure 12, now for P = 0.7. 
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Figure 14: (left): The residual uuf — Uiras of all points within 60° of the apex 
and antiapex of (l,b) = (0,90) and (300,0), for j3 — 0.4. (right) The same residual 
for the case (3 = 0.7. Note the linearly growing dipole residual toward (300,0). 
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Figure 15: The dipole field derived from the Mark III galaxies. The ITF dipole 
is plotted as the open symbols, while the IRAS dipole for /3 = 0.6 is plotted as 
the closed symbols. The X,Y,Z panels show the three Galactic components of 
the dipole and the top panel is their quadrature sum. The solid line is the full 
IRAS dipole field calculated from Equation 4. The Mark III galaxies are binned 
in redshift shells of 1000 kms~^ thickness. Note that the Ji?v4 5 inferred dipole for 
the Mark III symbols tracks the solid line, while the X and Y components of the 
ITF dipole are quite deviant from the line. 
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Figure 16: (top:) A scatter diagram of a**-^ versus the true coefficients a*^"'^ for 
the mock catalog. The ITF errors are shown, (middle:) Scatter of o:**'^ versus 
Q,iras f^j, i-j^g mock catalog. The increased scatter is induced by the errors in a*^"*, 
not shown because of the covariance. (bottom:) The scatter of a'*-^ versus a*^"^ 
for the real Mark III data (/3 = 0.6). Note the qualitative agreement of the mode 
amplitudes but the further increase in scatter. 
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